NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 





TECHNICAL NOTE 2378 


AUTOMATIC CONTROL SYSTEMS SATISFYING CERTAIN 

GENERAL CRITERIONS ON TRANSIENT BEHAVIOR 

By Aaron S. Boksenbom eind Richard Hood 

Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 





2093 


TECH UBRARY KAFB, NM 


] 



□□bS717 


NAHOHAL advisory COMMia?TEE FOR AERONAUTICS 


TECHNICAL NOTE 2378 


AUTOMATIC CONTROL SYSTMS SATISFYING CERTAIN 
GENERAL CRITERIONS ON TRANSIENT BEHAVIOR 
By Aaron S. BoksenLom and Richard Hood 


SUMMARY 

An analytic method for the design of automatic controls has 
been developed that starts from certain arbitrary criterions on 
the behavior of the contixiUed system and gives those physically 
realizable equations that the control system can follow in order 
to realize this behavior. The criterions used are developed in 
the form of certain time integrals. , 

General results are shown for systems of second order and of 
any number of degrees of freedom. Detailed examples for several 
cases in the control of a turbojet engine are presented. 


INTRODUCTION 

In the past several years, there has been increasing develop- 
ment and interest in automatic control; in the fields of gun 
direction, guided-missile control, and control of gas-turbine 
engines, for instance, where very refined and accurate controlled 
behavior is required, need still exists for further development of 
the methods of controls analysis and design. 

Recent developments in this field have been mainly concerned 
with the problem of control analysis both in the realm of linear 
systems (reference l) and in the realm of nonlinear systems (ref- 
erence 2) . These analytical works answer the following questions : 
How iri.ll a given system behave or how is its behavior affected 
when certain constants in the system are changed? 

Another problem of equal and, in some cases, greater import- 
ance is that of control synthesis. Work on this problem seeks the 
answer to the following question: Given certain criterions con- 

cerning the behavior of a controlled system, how should the system 
be designed? The answer to this question should give all aspects 
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of the system; for instance, whether the system shoTild he linear, 
"vjhat the general configuration should he, and what the precise 
values of all the constants should he. 

This synthesis problem has hardly heen broached in literattire . 
The use of analysis as a design procedure offers a partial solu- 
tion to this problem in that the analysis of a large number of cases 
may reveal, coincidently, a case that satisfies the desired cri- 
terions of controlled behavior. Such a method is, at best, long 
and tedious and almost always would result in conipromises because 
the systems chosen to be analyzed would probably be such that they 
could never satisfy all the desired criterions. 

A method for designing a linear system to satisfy certain 
special criterions when operating on a random input is developed 
in reference 3. This method is applicable aa an addition to a con- 
trol system xdienever random external disturbances are involved. 

A partial solution to the synthesis problem is developed in refer- 
ence 4 in satisfying the criterions of noninteraction for systems 
with many degrees of freedom. 

An analysis made at the HACA Lewis laboratory and presented 
herein develops a rational method of control synthesis that starts 
from any arbitrary but physically realizable criterions and results 
in the equations for the best system that satisfies these criter- 
ions. As is shown, the nature of the criterions, in generail, 
reqiiires minimizing certain time integrals by using the calculus 
of variations and the methods developed are an ^plication of the 
calculus of variations to the problen of control synthesis . 

A careful scrutiny is first made of the whole problem, fol- 
lowed by a development of general resiilts. These general resiilts 
are then applied in exanples to the design of turbo jet -engine con- 
trol systems. In general, the methods tised vary according to the 
order of the differential equation describing the plant, the proc- 
ess, or the system being controlled and according to the number of 
degrees of freedom or independent variables. Detailed analyses 
are presented for application to a system of first order and of 
one degree of freedom. General resxilts for systems of second 
order nrid any number of degrees of freedom are developed in the 
appendix. 
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SYMBOLS 

following symbols are used in this rejiort: 

gas-turbine-engine characteristic constants 

constant 

function of A 

functions of y and y 

indicates partial differentiation with respect to y 
indicates F [y(-ti), y(t]^^ 
arbitrary function 

. . weighting functions used for gas-turbine control 
criterions 

functions of y 

function used for gas-turbine-engine surge criterion 
function used for gas -tiirbine -engine blow-out criterion 
controller gain 

temperature limits for gas-tTorbine engine 

deviations of gas -turbine -engine speed, compressor- 
discharge pressure, characteristic tenperature, 
and fuel flow, respectively, from veilues at some 
common equilibrium condition 

actual engine acceleration 

exponent 

time 

time at end of transient 

variation in time at end of transient 

independent variables, functions of time 
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€ 

a 


T 


small, number 
arlDitrary constants 

gas-turtine-engine time constant for response to 
tenperature 

transient time constant of controlled system 


Subscripts : 

0 initial condition of variable 

f final value of variable 

s setting or desired value of variable 

Superscripts ; 

* indicates case different from optimum 

The dot indicates differentiation with respect to time. 

The prime indicates differentiation with respect to the argu- 
ment shown. 


ANALYSIS 

Survey of Problem 

Control problem. - An important aspect of the control -synthesis 
problem is a clear definition of the criterions of desired con- 
trolled behavior. If a variable y is to be controlled, a reason- 
able criterion is that the time integral of some function of y 
is to be a minimum or a constant; that is. 



f(y) dt = constant or minimum 


( 1 ) 
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or 


ntj 


\J 


(y-yg)^ 


dt = constant or minimum 


( 2 ) 


Equation (2), for instance, weigSits the error in y eis the square 
and according to the time duration of that error. Another type of 
criterion may he that Tdilch requires a certain time duration to he 
a minimum or a constant; that is. 




U 


0 


dt = constant or minimum 


( 3 ) 


The use of a single criterion, such as equation (l), vill \isu- 
ally yield f(y) = constant. This result is reasonable because 
f (y) can tistxally he made identically a constant if no additional 
criterions are inposed on other variables in the system. Usually, 
certain limiting conditions exist, however, on other variables in 
the system and these conditions must be included in the original 
criterions . 


Thus, a possible criterion could be written as follows: 


for 


n ‘^1 


(y-yg) = minimum 


O'^l 




f(z) dt = constant 


} ( 4 ) 




If, for instance, y = engine speed and z - characteristic tem- 
perature of a gas-tiorbine engine, the criterion of equation (4) 
states that it is desired to design a control system such that, 
for a particular value of a temperature integral, the integral of - 
the speed-errof squared is a minimum. This criterion may be used 
if, for instance, it is known that an over-temperature condition 
can be tolerated for a certsiin period of time and it is desired to 
keep the average speed error at a minimum during a transient. 
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The general theory vill show that as many crlterions eis 
desired of the type «howa in eguations (l) to (4) can he included 
together and a control system can he derived that automatically 
satisfies all these criterions simultaneously. 

Another aspect of the control criterions is the end conditions 
of the Integrals of equations (l) to (4) . The time interval for 
■which these integrals are to he a minimum or a constant must he 
chosen. A reasonable time interval is any duration during which 
essential ejctemal dlsturhances are constant and during •which the 
system to he controlled moves from one essential level of opera- 
tion to another. The essential external dlsturhances are those 
that cannot he immediately corrected hy the control system. If an 
essential external disturhance were allowed •within the time inter- 
val of the criterions, no physically realizable system could he 
expected to anticipate this diaturhance so as -to behave properly 
before this disturhance occurs. An essential level of operation 
is any specific condition of only those variables that must he 
continuous. It -vri-ll he shown that the essentieO. level of opera- 
tion appears as a natural houndaary condition for the type of cri- 
terion used. In the case of a turbojet engine the transient 
heha^vlor of -which can he described hy a first-order differential 
equation, -fche engine speed determines the level of operation. If 
a lag exists in the fuel system or he-tween tempera-ttire and engine 
speed, then both engine speed and acceleration are required -to des- 
cribe the essential operating level of the engine. 

Analytic problem. - The con-brol system resulting from any 
design method must he physically realizable. Th^e are two aspects 
to -this problem. First, it is possible to set do^wn criterions that 
are not realizable -with any system or are incompatible -with each 
other. If such criterions are •used, the unrealizability -will 
appear either as a requirement on the control to look ahead into 
the future or as an inability to satisfy the boundary conditions 
of some differenticLl equation. In most cases, a clear iinderstand- 
ing of the criterions used and of the system to be controlled -will 
indicate inconpatibilities of this sort. 

The second aspect of physical realizability is pxirely mathe- 
matical. It is desired to derive a description (a differential 
equation) of the control or the controlled system that satisfies 
the criterions of con-fcrol and all the necessary boundaiy ^conditions 
that arise in the derivation of this equation. Although the ma-the- 
matical solution of the problem may be any derivative or integral 
of -this differential equation, -fche physical solution of the prob- 
lem requires the differential equation that itself satisfies the 
boundary conditions and for which no -undetermined constants of 
integration exist. Thus, such forms as 
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and 


y = Cx 


y = Cx 



(5) 


are not necessarily interchangeable as descriptions of some part of 
a controlled system because the forms differ by an undetermined 
constant of integration. For stable linear systems, the effect of 
this constant becomes vanishingly small; for the general nonlinear 
systems presented herein, however, this constant must be considered. 


Stability problem. - The requirement of stability is a special 
criterion that does not enter into the main body of the methods of 
this report. It may enter in the final steps of the method where ' 
the final differential equation describing the controlled system 
may be the integral of a higher-order differential equation that 
satisfies the necessary boundary condition for stability. In addi- 
tion, it is always necessary to add to the controlled system a 
stability device that does not affect the behavior of the system 
as far as satisfying the other criterions is concerned. This 
device can be described as follows: 


when 

y = yg 

then 

y = 0 

or, for a second-order system, 
when 

and 
then 




( 6 ) 


( 7 ) 


General Theory and Results 

It has been shown that the criterions for control can be 
developed in the following forms: 
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riti 


U 


f(y) at 


0 




U 


(y-yfi) dt 


f^(z) dt 





) ( 8 ) 


y 


and so fortli. If, for such a 3J.st of criterions, one of the inte- 
grals is to he a miiilmum under the condition that the other integrals 
are to he constant, it is sufficient (reference 5) to make 



The A’s are arhitrary constants that enter into the control sys- 
tem as the adjustable parameters and are precisely deterMned hy 
the choice of values that the constant integrals are to have. 

The technique of the A multipliers is -widely used for prob- 
lems of this type -where one condition is to he a minimum under 
other restrictive conditions. Indeed, the conditions need not he 
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in integral form and any functional or differential relation among 
variables can be handled in a similar manner (reference 6) . 

Equation (lO) can be made very general "when all possible 
restrictive conditions are included. In the final eq.uations, -which 
are derived later, if any one criterion is not to be used, then the 
corresponding X — ►O. If any of the criterions are -fco be zero (as 
for the case of a variable having an absoltrte limit), then the cor- 
responding "K — » 

If, for this development, the system to be controlled is of 
first order and of one degree of freedom (has one independent vari- 
able) , then the variables y and z are such that z = z(y,y) . 
Equation (lO) can then be -written, in general, as 


F(y,y) dt = minimum 


In equation (ll) , F is a continuous function of y and y, and 
y is a continuous function of time. 

The calcultis of -wnriations (reference 5) is used to determine 
y as a function of time such that the integral of equation (ll) 
is a minimum; that is, if the solid curve of figure 1 makes the 
Integral a minimum, any other curve (such as the dashed one) -will 
make the integral equal to or greater than the first integral. If 
•the t-wo curves are very close, the condition that the t-wo inte- 
grals are equal will make the integral of equation (ll) stationary 
(maximum, minimum, or Inflection point). 

The problQtt is similar to that of finding the mH.-iriTmnn or mini- 
mum point on a curve by setting the first derivative equal to zero, 
^^ether a maximum or minimum point exists- is decided by the second 
derivative at that point. In the variational problem, pro-vLng a 
true minimum involves taking the "second variation." As in -fche 
problem of finding a minimum on a curve, the second variation proves, 
at most, a local minimum and not an absolute minimum. For the 
specific examples discussed herein, an absolute mi-niimm is proved 
by another method. In many cases, the physical meaning of the 
equations -will indicate the existence of a -unique mi-nlTmnn that 
obviates going further than the "first variation." 
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On figure 1, the curve for minimizing the integral of equa- 
tion (u) is the solid curve having the value y(t) at any time. 
Any other curve differing from it hy a small amount is shown as the 
dashed curve having the value y(t) + € 5y(t) at any time, where 
e is a small nuniber and 5y(t) is an arbitrary continvious func- 
tion of time. The' condition for the integral of equation (ll) to 
be stationary is 


ptj^+e StjL 



F(y+£ 5y, y+e Sy) dt = 0 


( 12 ) 


at 


e= 0 


■where — [5y(t)3 = By(t) (reference 5) . The time duration for 
dt 

the integral of equation (12) is such as to start at some definite 
time (t=0) but not to end at a definite time but rather along 
some curve y = f(t) (fig. l), in order to allow the proper bound- 
ary conditions of moving from one essential level of operation to 
another, as previously discussed. 


Performing 



Integrating the 



the operation indicated by equation (12) leads to 
^1 

Fy 6y dt + I Fy 6y dt + F(ti) bt^ = 0 (13) 



second term of equation (13) by parts gives 


It-, 


- ^ (Fp] 87 dt + Fy By 


+ F(ti) btj^ = 0 (14) 


Because 6y is an arbitrary function, the integral and the 
boundary-condition terms must vanish separately. The geometry at 
the end condition (t=tq, fig. l) gives 
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Sy(ti) = |f‘(ti) - y(ti)] 8ti (15) 

as € —*■ 0 . The two conditions that follow from equation ( 14) then 
become 

J ’^’^1 


and 


5ti |F(ti) + Fy(ti) [f(ti) - y(ti) I - Fy(0) 5y(o) = 0 


(17) 


The time interval d\iring "vdiich the criterion of equation (u) 
is to hold is considered as that daring which the system moves from 
one essential operating level to anotherj in this case, from one 
definite value of y to another definite value of y. Thus, 


5y(0) = 0 


and 


f’(ti) = 0 


(18) 


Eq.uation (l6) is satisfied only if the integrand is zero, and, 
because bt^^ / 0 , the two conditions of equations (16) and (17) 
become 



and 


F(t^) = Fy(t^) y(t^) 


(19a) 


(19b) 


Equation (l9a) need not hold at t = 0 because &y(o) = 0. The 
only condition that need hold at t = 0 is that Fy(0) is finite, 

tdiich will be true if y is continuous at t = 0. At the start 
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of a new •transient, y, F, F , and F* may te discontinuous, 

•wiiereas at other points (O < t < t^^) F^ will he continuous 
because of equation (I9a) . 

Equation (l9a) is the differential equation for the y(t) 
that satisfies the oidginal criterion of equation (ll) . The physi- 
cal answer -to the problem is the first integral of equation (l9a), 
which satisfies the boundary-condition' equation (l9b) . This solu- 
tion is 

F(y^y) = y Fy(y>y) (20) 

and wheneTer y, Fy, and so forth, are continuous. 



by differentiating equation (20) with respect to time. Thus, 
either y = 0 (which is only true during stability) or equa- 
tion (l9a) is satisfied. 

Thus, equation (20) is the description of that physically 
realizable system the behavior of which will automatically and 
simultaneously satisfy those criterions included in the fimction F 
during that time interval for which the external disturbances are 
constant and during which the system goes from one operating level 
to any other operating level, A stability device must be added to 
■the systemj the description of such an ideal device is as foUovre: 

when 

y = ys 

then 

y = 0 

An additional condition must be met if y is discontinuous in the 
interval 0 < t < t^^. In this case, Fy must be contimwus during 

the discontinuity in yj Fy, however, -will usually be discontinuotis 
■vdien y is discontinuous. This discontinuity in Fy usually 

means that some essential external disturbance has entered the sys- 
tem and the point of discontinuity imist be the start of a new time 
interval. 
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Application to Turbo Jet -Engine Controls 


In the usual case of designing turbojet-engine controls, the 
engine speed that sets the essential operating level of the engine 
and, in the main, other pertinent characteristics such as thrust, 
is to be set or controlled. Limiting conditions of the engine are 
those of overspeed, overtenperature, compressor surge, nrifl rich 
bximer blow-out. The following criterions on the behavior of this 
engine are typical: 

^t^ 


f^(N-Ng) dtj for speed control 


f«(R) dtj for speed overshoot 


fgCT) dtj for tenperature overshoot and 
undershoot 


> (23) 


^ti 


- g(N)] dtj for compressor surge 


fgjp - h(N)^ dtj for blow-out 


dtj for rise time 




Figures 2(a) to 2(e) show the nature of the fimctions appearing in 
equation (23). The variable involved is essentially weighted by 
the fimction shown and integrated with respect to time. The quan- 
tity P - g(U) is the amoxint the compressor-discharge pressure 
exceeds the safe pressure for surge and g(E) (shown in fig. 2(f)) 
is the compressor-discharge pressure for each engine speed at a 
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safe value ‘below surge. Eicli ‘burner ‘blow-out can ‘be handled in a 
similar manner (h(N) is shown in fig. 2(g)). The rise time is 
the total time for the system to move from one essential operating 
level to the other. 

The linearized engine characteristics can he expressed, hy 
assuming first-order behavior, as follows: 

T = aN + aoK (24a) 

P = hU + cT (24h) 


Thus, the integrals of equation (23) are all of the form 

J ■^tl 

f(N,H) dt (25) 

a 

where f is a continuous fxmction of H and K, and H (barring 
impulsive tenperatrires and the like) is a continuous function of 
time. 


Speed control; case A. - If only the enror in speed control is 
considered important, the criterion becomes 


fj^(N-Ng) dt = minimum (26) 

where F = fi(N-Ng) . Equation (20) becomes 

fi(?r-Ng) = 0 

and from the nature of f^^ (fig. 2(a)), 

N = N 
s 

This result means that, in the absence of other criterions on the 
engine behavior, this speed control should keep speed error identi- 
cally zero, which is physically realizable only in the sense of 
allowing infinite tenperatures and the like to keep the speed error 
identiceilly zero. This result, however, is inconsistent with the 
previovis development of equation (20) .in that N is now a 
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discontinuous function of time and the time intervgil of the inte- 
gral of equation (26) is zero. This instance is actually a trivial 
case of the general problem. The result {equation (27)) does indi- 
cate that a criterion like that of equation (26) must be acconpanied 
by an additional criterion (tenperatvure overshoot, for instance) 
to give a physically realizable system. 


Speed control uith temperature-limiting criterions; case B. - 
If the error in speed control and the overshoots and undershoots 
in teniperature are to be considered as the primary criterions of 
control, then from equations (lO) and (U), 


pti 


vjo 




dt = minimum 


(28) 


where F = f , (N-K ) + A f_(T). Equation (20) then becomes 
J. s o 

fl(N-Ng) + A f3(T) = A aoK f3>(T) (29) 

and the ideal stability device is such that 


when 


N = Ng 


then 



(30) 


Equations (29) and (30) describe the coc 5 )lete control system. 


In figure 2(c), it is convenient to let fs(T) = (T-L 2 ) for 

T > Lg and f 3 (T) = (T-L 

should be >1, because, when n s_l, T may be infinite and of such 
nature as to make N discontinuous and physically unreal even 


) for T < Lj^. In general, the power n 




though ^ 

and fi(N-Ng) = 


is finite. 

(N-Ng)^. 


In the examples of this report. 


n = 2 
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A variety of methods of setting up the control system to 
ize equation (29) exist. When the preceding egressions for 
and. f 3 are used., equation (29) can he put in the convenient 

— + (L-aN)^ = 

A 

•where, for eicceleration, •when K < Hg, then 

N > 0 and L = Lg [ 

and, for deceleration, •t.dien N > Hg, then 

H < 0 and L = In, 

J 

A schematic diagram of the system is sho^wn.in figure 3, where 
equation (3l) is considered to give a desired N for any value 
of N, Ng, and L. Consistent -with equation (3l), a rigtt -triangle 
construction is -used, -bo give a desired U. The ac^tual N can he 
made very close -fco N hy -using a high-gain proportional controller, 
as sho-wn. Pro-vision must he made to change the sign of N and 
the value of L -when K-Ng changes sign. In addition, the sta- 
hili-fcy condition requires a pro-vision for making if = 0 whenever 
speed error is very small or zero. 

The -use of a high-gain .proportional controller, -which follo-ws 
from the reqiiirement that N may he discontinuous, means that the 
fuel-flo-w rate required may he infinite if lags exist in the fuel 
system or in the feed-hacks. But, as no criterion heis heen set on 
fuel-flow rate, this requirement does not -violate the original 
criterions. If necessary, ho-wever, a criterion on fuel-flow rate 
may then he added -bo equation (28). Even though a criterion on 
tenperature is being satisfied, figure 3 does not use any direct 
meastirement of tempera-bure . ActueOly, the equation for tempera- 
•bure (equation (24a)) is used as an indirect indication of 
teDpera-bur e . 

The control system of figure 3 has one adj^lstahle parameter A. 
For any value of this system -will, for the value of integral 
tenpera-bure-overshoot squared obtained, give the mln-fTmjm value of 
integral speed-error squared. The value of A determines the 
ac-bual value of the integral tenpera-bure-overshoot squared. 


real- 


form 


(31) 
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. The integral ten^eratur e-overshoot squared as a function of A 
is shown in figure 4 for the special case where aJJg = Lj that is, 
acceleration or deceleration to the speed that corresponds to lim- 
iting teniperature . In this case, the system of figure 3 becomes a 
siu^tle first-order lag system and equation (3l) becomes 

E(L-sN) = aoH (32) 


•vdiere 



In figure 4, the Integral speed-error squared, the TOaTH-rmnn 
teniperature, and the time constant for this transient are also 
« shown as functions of A. A curve showing the miu-tTmni rise time 
for the corresponding temperature integral is included for compari- 
son that will be discussed later. The equations for the curves of 
figure 4 are as follows: 



(E-1)^ 

2E 

- 2l. 

~ 2E 


1 


(time constant) ^ ~ ^ 


> (33) 


The left side of these equations have been put in dimensionless 
form. The maximum tenperature T^^ny occurs at the beginning of 

the transient. The time constant t is that for the controlled 
system and is shown compared with the engine -time constant a. 

From figure 4, the value of A is chosen as a conpromise 
between the various quantities of equation (33). For E = 1 (A=o>), 
the temperature does not overshoot, the speed integral is 0.5, and 
T = 0 . As E increases ( A decreasing) , the temperature integral 


0 
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and the mflTrtTmnn tenrperattire increeise, -whereas the speed integreO. 

pnrl the time constant decrease. A value of E = ^\f2 (a^A = l) 
appears "to he a good conpromise anfl is -used for, the subsequent 
discussion. 

The heha-vlor of the system of figure 3 can be best seen by 
drawing the -tarajectories in the phase plane sho-wn in figure 5 , 

Inhere aoE is plotted against aU for lines of constant aNg 
aiccording -bo equation (3l) . Lines of constant tenpera-ture are 
45° parallel lines in this plot and the lines of T = L are shown. 
Each trajectory intersects, and is tangent "to, the line T = L at 
ET = Hg. Figure 5 completely describes the transient behavior of 
the system. For any starting point anywhere on the plane (for 
instance, point A), the system -will au-tomatically move the oper- 
ating point to that trajectory corresponding to the Ng that 
exists (point B) and then along this trajectory to the point T = L, 
(N = Ng, point C), and finally the stability condition will enter 
to move the operating point along the solid vertical line to 
N = 0 (point D) . 

The time sense for these transients is obtained by sol-vlng 
■the differential equation (equation (3l)) for -the speed and the 
temperature time responses. 53ie equations for these solutions are 
as. follows; 



(34) 


These transients are shown in the general case in figure 6 for a 
step increase in Ng. Maxiinum tenpera-ture occurs at t = 0 and 
the tenpera-ture overshoot decreases as N increases and when 
N = Ng, T = L and T = 0. The stability condition then causes T 
to drop to its equilibrium value at this point. The time scale 
sho-wn in figure 6 corresponds to the ^ecific relative values of 
•the ordinates shewn. 

M-t-nirmnn rise time and tenpera-bvure-H-Tniting criterions; case C. - 
In order to obtain a Tnlu-tminn rise time for a constant tenpera-ture 
integral, the requirement is that 


0 
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[i + x fjd)] 


dt = minlTmuti 


•where F = 1 + 7vf3(T). Equation (20) hecomes 


1 + ^ fgCT) =AaaN fj'CT) 


(35) 


(36) 


and the ideal stahility de'vlce is such that 


when 


then 



(37) 


•where Ng is the desired stable N. Equations (36) and (37) des- 
cribe the complete control system. Eq.uation (36), and therefore 
the control system proper, does not include Ng and, in this case, 
the stability de^vlce is independent of the control system proper. 


\Jhen the same f3(T) is used as -was oised in the previous 
case, equation (36) can be put in the convenient form 


^ + (L-aK)^ = a^a^ 

for acceleration, -when If < Wg , then 
N > 0 and L = Lg 
)r deceleration, when N > Eg, then 




(38) 




A schematic diagram of this system satisfying equation (38) is 
sho^wn in figure 7. This control system is the same as the pre^vious 


system (fig. 3) except that in figure 7 


replaces 


H-Ns 

aJ a 


on 
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one leg of the right-triangle construction and the stability con- 
dition must be iiogposed outside the control system, as Hg does 
not enter into the criterion of control (equation (35)). 

The system of figure 7 has one adjustable parameter 7^ -which, 
as before, sets the precise values of the integrals entering into 
the criterion (equation (35)), as -well as all other beha-vior 
characteristics . 

The tenpera-ture integral is sho-wn in figure 8 as a function 
of A for the specieil case -where aNg = L; that is, for accelera- 
tion or deceleration to the speed that corresponds to li m iting tem- 
perature. In addition, the rise time, the ma.yl.mum tempera-ture, 
and the integral speed-error squared are sho-wn for this transient. 
A cmrve pho-wing the Tn-tTiimiim speed integral for the corresponding 
tenperatu.'e integral (from case B) is also sho-wn and -will be’*dis- 
cussed later. The equations for the curves of figure 8 are ejs 



(39) 

Again ‘A is to be chosen as a conpromise be-tween these various 

1 

quantities. In this case, units of temperature 

and V^(L-aNo) is the pertinent dimensionless parameter. For 
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a fixed "Kt the initial speed Nq ¥111 therefore change the time 
integrals and an increased L-alTQ vill decrease the tenperatnre 

integral and the Tna-idmiTm temperature and increase the rise time 
and the speed integral. 

The dynamic hehavior of the system of figure 7 is^shoTm in the 
phase plane in figure 9. This figure is a plot of aoN against aH, 
according to equation (38), for various values of X . Lines of 
constant tenperature are 45° parallel lines and the lines T = L 
are shovn. For any fixed A, only two trajectories are followed: 
one for accelerating and one for decelerating. From any starting 
point on the plane, the system will automatically move to that 
trajectory corresponding to acceleration or deceleration and will 
move along this trajectqry until the st ability condition enters 
(at K = Hg) to make N = 0. 

The dependence of the time integrals on Nq may require X 
to vary with Nq* In figure 9, Ho*s corresponding to each 7^ are 
shorn such as to keep A(L-alfo)^ = 8. The value of this parameter 
was taken so as to have the temperature integral in this case 
equal to that of the previous case for pxirposes of conparison. 

The time sense for these transients is obtained by solving 
the differential equation (38), for idle U-t and T-t transients. 
The equations for these solutions are as follows: 



These transients are shown for the general case in figure 10 for a 
step increase in Ng. The tenperature will jimp to some value 
above L, but in this case, unlike the previous case, the tenpera- 
ture continues to increase as the speed Increases. Maximum tem- 
perature occurs at the end of the transient. IJhenever N = Ug, 
two such conditions are shown, the stability condition takes over 
and T is reduced to its equilibrium value. The time scale shown 
on figure 10 corresponds to the specific relative values of the 
ordinates as shown. 

Gonparlson of cases B and C . - In case B, a system was devised 
that, for a constant value of the integral temperature-overshoot 
squared, gives a minimum integral speed-error squared. In case C, a 
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system was devised that, for a constant temperature integral, 
gives a mini mum rise time. For case B, the rise time is not a min- 
imum hut can he congiared with the minimum rise time of case C . 

For the specisil transient of accelerating to T = L, the system of 
case B reduced to a first-order system the time constant of which 
is shoim in figure 4. Because five time constants are considered 
as the rise time of an exponential, the time constant of ceise B 
should he conpared with one-fifth the miniirnmn rise time of case C . 
In figure 4, for corresponding values of the tenperature integral, 
one-fifth the minimum rise time from case C is plotted. Figure 4 
now shows that the rise time for case B is about twice the ml.nlTnum 
possible rise time. 


At corresponding values of the tenperature Integral, 


(W-Nf ) dt for case C can he compared with the minimum 


possible value of this quantity from case B. Ifhe mi nimum veilue 

lti 

2 

of / (N-N^) dt is plotted in figure 8. It is seen that, for 



the same temperat'ure integral, I (W-N^) dt when minimum 


rise time is obtained is about 120 percent of the minimum possible 


value of 


(N-N„) dt . 


Proof of absolute minimum. - As previously noted, the minimi- 
zation of equation (H) involves not only the condition that the 
first variation be stationary, idiich leads to equation (20), but 
also that the second variation be positive. This second condition, 
however, would prove, at most, a local minimum for the condition 
of equation (U) . A special proof of an absolute minimum is shown 
616 follows: 


If equation (ll) is written in the following form: 


F(y,y) ^ = minimum 
y 



( 41 ) 
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and y is considered as a function of y for the minimum condi- 
tion and ^ is considered as a function of y for any other pos- 
sible cane, then the condition for an absolute minimum is that 



>rtiere dy/y* is a positive differential. If the function F is 
considered in a general form^ quadratic in y, as foUovs: 


F(y,y) = G(y) + y G^Cy) + y^ GgCy) (44) 


then equation (43) becomes 


pyf 

r* *1 

1 (y*-y) 

^ G 2 (y) - ^ 

Jyo 

y _ 


7 * 


>0 


(45) 


The equation of the control system (equation (20)) becomes, 
for the form of F assumed in equation (44), 


.2 


y Gg(y) = G(y) 


(46) 


Using this e3q)ression for G(y) puts equation (45) in the form 

lyf 



•-2 ay 


&2(y) (y*-y) p - ° 


(47) 


p p 

For cases B and C, G2 = Aa^a“^ -where 0. Thus an absolute mini- 

mum is proved for these t-wo cases. 
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Degree of ml.nlTiium, casQ B. - The left side of ec[uation (47) 
can he used as a measure of the deviation from optimum conditions 
uhen the methods of this report are not used. If the two cases 
(W for .the optimum case and for any other case) are compared 

at the same value of the tenrperature integreQ., the left side of 
equation (47) is the difference between the integral speed-error 
squared for any case and the Tni.ni-mnm possible value of this quantity. 
The ratio of this deviation to the mlnlimmi value becomes 


Fractional increase 
in speed integral 



(48) 


for the transient of acceleration or deceleration to limiting tem- 
perature. The coefficient of the integral of equation (48) for 

the value of 'K previously chosen (a^A= l) is 2\[2. 

Case B can also be considered as giving a mi-nlTmrm value of 
the temperature integr^ for any definite value of the speed inte- 
gral. If two cases (N for the optimum and for any other 

case) are compared at the same value of the ^eed Integral, the 
left side of equation (45) is proportional to the difference 
between the integreil teniperature-ovei'shoot squared for any case 
and the minimum possible veOue of this quantity. This deviation 
can be written .as 

p^f 

Fractional increase in _ 2E | /o*-T \ dy 

temperature integral (E-l)^JyQ ^f'^q/ ^ 


for the transient of acceleration or deceleration to limiting tern- ^ 
perature. The coefficient of the integral of equation (49) for 

the value of previously chosen (a^A = i) is 16.4. 

Degree of mi.nitmim, case C. - If the two cases N and K* are 
con 5 >ared at the same value of the temperature integral, the left 
side of equation (47) is the difference between the rise time for 
any case and its minimfum possible value. This deviation can be 
Witten as 
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Fractional increase 
in rise time 


■K (L-aHp)^ 

sinli"^V>v(L-aNo) 



(50) 

for the transient of acceleration or deceleration to limiting tem- 
perature. The coefficient of the integral of equation (50) for 
the value of /s/^(L-aUQ) = 2*</2 previously chosen is 4.54. 


Case C can also he considered as giving a mininunn value of 
the teD5>erature integral for any definite value of rise time. If 
the two cases H and are now con^iared at the same value of 

rise time, the left aide of equation (47) is proportional to the 
diJPf erence between the integral temperature-overshoot squared for 
any case and the ml-nimum possible value of this quantity. This 
deviation can be written as 


Fractional increeise in 
temperature integral 



^/>v(L-aIIo)2 + 1- V^L-aHo) 


(51) 


for the transient of acceleration or deceleration to limiting tem- 
peratTire. The coefficient of the integral of equation (51) is 16.4 
for the value of /\/T(L-aITo) previously chosen. 


GENERAL SUMMARY 

When the criterions on the behavior of a controlled system can 
be expressed in certain general forms, 6is follows: 
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r»tn 


U 


f(y) dt 


0 


(y-ys) dt 


dt 


dt 




) ( 8 ) 


J 


where the time interval is taken as any duration d\iring which 
essential external disturhances are constant and during which the 
system moves from one essential operating point to another, the 
optimum control can he considered as one that minimizes one of the 
integrals of equations (8) while maintaining the other Integrals 
at prescribed values. The analytical problem, according to the 
calculus of variations, reduces to the following equation: 


r*i 


UO 


f(y) dt + I (y-Yg) dt + ^2 I dt + X3 I dt =■ minimum 


( 9 ) 


For general first-order systems, where z = z(y,y), equation (9) 
reduces to 


F(y,y) dt = minimum 


( 11 ) 


where F = f(y) + Tv^Cy-yg) ^3* equation neces- 

sary for the control system to satisfy equation (h) and all the 
boundary conditions becomes 
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This equation should he followed hy the control system proper. In 
addition, a stability device must he added to the system the ideal- 
ized characteristics of which wotild make y = 0 when 7 = y^. 

The arbitrary multipliers \ are then found hy evaluating 
the integral criterions involved in F- The transient behavior of 
the system derived is found hy Solving the differential equa- 
tion (20). The degree of the miniTnum or the amount suffered when 
any other control system is used was evaluated for the special 
cases considered. A summary of these developments follows for a 
special form of the F function, quadratic in y, -tdiere 

.2 

F(y^y) = G(y) + y G^^Cy) + y G2(y) (44) 


Control-system eqtiation. - For F in the form of equa- 
tion (44), the control-system equation (20) becomes 


.2 G(y) 

&2(y) 


(52) 


The function G^^(y) does not affect the control system and the 

control-system equation gives an explicit expression for y els a 
function of y. 

Evaluation of integrals. - If an integral of a func- 
tion H(y,y) is any one of the criterions to be considered, it can 
be evaluated els follows: 



Thus the integrals can be evaluated ifithout solving the differen- 
tial equation (52) . 

Transient behavior. - The differential equation (52) can be 
easily solved as follows: 



( 54 ) 
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Degree of Tnl-nimum. - Equation (47) 'vreis derived to prove the 
absolute mln-iTTiinTi for the special cases considered and can he used 
to evaluate the degree of the Mnimum found. Thus, 




rt-, 


H(y , dt* - 
0 UO 


H(y,y) dt = 


T 


pti 


■v^ere 


UO 


GgCy) [^(y) y(y)J 

(55) 


H(y,y) dt is to he a Tni.nimnn or maxiinum. 


If ^ ^ ^0 ^ ^ ^f ' then an absolute TtiiniTnuni is 

obtained; if G 2 (y) < 0 for y^ < y < y^, an absolute maximum is 

obtained. Equation (55) also Indicates that the degree ot minimi am 
(or maxiraum) varies -with the magnitude of 02(y) * ^ ^ 

denominator indicates that for any deviation j^-yj in y, it is 
better to err on the side of a larger |* 


SUMMARY OF RESDITS 

A rational analytic method for the design of automatic con- 
trol systems has been derived. Criterions on the behavior of the 
controlled system were developed in the form of certain time inte- 
grals. When any of these arbitrary but physically realizable cri- 
terions were used as a starting point, those equations that the 
control system must follow were derived. The criterions developed 
required the minimization of certain time integrals using the cal- 
culus of variations. The method gave not only a description of 
the behavior of the controlled system but also gave those physically 
realizable equations that the control system can follow in order 
to realize this behavior. General results were shown for systems 
of second order and of any number of degrees of freedom. 


Leirf-s Fli^t Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, October 11, 1950. 
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APPENDIX 

SECOND-ORDER SYSTEMS OF SEVERAL DEGREES OF FREEDOM 


It is beyond tbe scope of this report to detail the cases of 
higher-order systems and those having several independent variables. 
The general equations for these cases -will be developed and it is 
ejtpected that subsequent reports ■will cover their applications in 
detail. 


For the case of a second-order system vith two degrees of 
freedom, equation (ll) becomes 


Oti 


Uo 


F(y>y^y^z^z/z) dt = minimum 


(Al) 


where y and z are independent functions of time. The condition 
to satisfy equation (Al) is 


pit^^+G&t]^ 


de 


U 


t=o 


F(y+e5y,y-fc£5y,y-K:5y,z+£5z,z+e5z,z+e6z) dt = 0 (A2) 


at 


c= 0 


The time duration of the integral of equation (Al) is that begin- 
ning at a definite time (t=0) but not ending at any definite time 
but rather along some curves: y = fi(t), y = f 2 ('t), z = g^(t), 

and z = g 2 (t) (fig. H) • The functions 6y and 5z are arbi- 
trary and independent functions of time. 


Performing the operation indicated by equation (A2) gives 



(Fy 6y + F^ 6y + F” 6y + F^ 6z + F* 5z + F” &z*) dt + F(tj^) 5t^ = 0 


(A3) 
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After integration "by parts. 


nifc]_ 



0 


ir 2 1 


(Fy 6y + Fy 5y - ^ Fy 5y)j^^ + FCt^) Bt^ + (A4) 


[f- 5z +F”5z - 


= 0 


As tefore, the integrands of the integrals and the houndary- 
condition terms must vanish separately. From the geometry at the 
end condition (t=t^, fig. H)j 


By(\) = 

6y(ti) = ^^1 

&z(ti) = [gi’(ti) - z(ti)] 5% 

5z(ti) = [g 2 ’(ti) - z(ti)] 8ti 

The three conditions from equation (A4) then become 




(A5) 


J 
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'N 


''y - K ^ <^y) ■ ° 


r - yTy ■'■ y ^ (^y) ■ y^y ■ ^2 ■'■ ^ ^ ' ’^ 2 ] 

^i'(ti)[j'y - ^ yy(%) + 

Si'(ti)jF5 - ^ + 

'Ify ■ « ^^y*]t=o ^ * 


t=tn 


>(A6) 


6y(0) 


5z(0)|f- - A (Fv)j^^ + 5z(0) F”(0) = 0 


y 


The first two of equations (A6) are the differential .equa- 
tions that satisfy the original criterions of equation (Al) . The 
physical solution to the prohlem is the pair of solutions of these 
equations that satisfy the houndary-condition equation of equa- 
tion (A6) . If the first of equations (A6) is multiplied hy y 
and the second hy z and the equations are added, an exact deriv- 
ative is formed the integral of -which is as foUo-ws; 


F - + y (FV) - ’ 3 ^” - ^* + z ^ (FV) 


dt 


dt 


zF^ = C 


(A7) 


A s-tudy of the houndary-condition equation of equation (A6) then 
sho-ws that the physically reasonable houndary conditions- shoTild he 
as follows; 



32 


NACA TN 2378 




If 

^y 

If 

^y 

If 

F* 

^ z 

If 

F" 

Then 


If F" / 0, then 6z(o) = 0 and g 2 '(t 2 ^) = 0. 


J 


(AB) 


JLXX u wxo. \Ai } y w ~ w & wuw -l. i, i ir*ju w« < 

prohlem of equation (Al) for the boundary conditions of eq.ua- 
tion (as) is as follows; 


F - :^y + y ^ (Fy) - ‘jFy - 2 ^ (^z) - ^ 0 




and either 


or 


p-o- - — (F-i;-) + ^ = 0 

y dt y ^^2 '• y' 


F _ A- (f*) + — (F'O = 0 
it '' 2^ ^ ^.^2 


> (A9) 


The meaning, of the boundary-condition eq.uation (A8) is to 
define the original criteriohs for that duration during -vdilch the 
system moves from one essential operating level to another. Thiis, 
if all conditions of equation (AB) must hold, the system goes from 
one definite y, y, z, and z to any other definite y, y, z, 
and z. The first differential equation of equation (A9) would be 
of third order and the second or third equations of (A9) woTild be 
of fourth order. If equations (A9) are integrated to obtain a 
pair of second-order differential equations having three constants 
of integration, the choice of these three constants can then deter- 
mine a desired end point; that is, the values of y, z, and z at 
some final y. 

The physical solution to the problem, then, is the pair of 
second-order differential equations that are solutions of equa- 
tion (A9) and the constants of integration of which are evaluated 
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so that the system goes through some desired end point. This end 
point must he such as to allow the possibility of stability. Such 
an end point may be written ais follows: 


when 


then 


y = yg 




y = 0 

Z = Zg 

z = 0 


/ (AlO) 


J 


which- gives three conditions for the e-valuation of the three con- 
stants of integration. A stabili-ty de-vice must still be added to 
the system so -that, at the point -when equation (AlO) holds. 


and 


• • ^ 

y = 0 


*z = 0 



(All) 


Equations (A9) are symmetric in the variables y and z, 
■which indicates the nature of the extension for irnsre independent 
variables. Thus, if a third independent function w exists, the 
original criterion -woiild be -written as 


F(y>y/y^Z;Z,z',w,-w,w) dt = minimum 


(A12) 


This condition is satisfied under boundary conditions similar to 
equation (AB), where two additional conditions are added on the 
variable w and the folio-wing equations describe the controlled 
system: 
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F - jF* + y — (F*0 - W** - zF* +. z — (F**) - zF** - 

y dt y y z at z z 


"hF ' + w — (F* - tF* ’ = 0 

and any tvo of the folio-wing three equations: 


S H = 0 

dt 




>(A13) 


J 


The three equations of equation (A13) can then he integrated to 
give three second-order differential equations where the five con- 
stants of integration are evaluated so that the system goes throu^ 
the desired value of y, z, z, w-, w for some final -value of y. 
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fl(N-Ns) 





(c) Temperature overshoot 
and undershoot. 




below surge. 



Figure 2. - Arbitrary weighting functions for various control crlterlons 
and pertinent engine characteristics Involved. 
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Pigiire 3. - Schematic diagram of control for turbojet engine for 
case of speed control with temperatxire-llmlting crlterlons 
(case B) . Stability device must be added to this system. 



Temperature integral = — — 

0 „2rN^_M 

^ — x--v^ -- 

1 r 

speed Integral = — r- 

^ o (Nq-Nb) 

V* 


1/5-iuln. rise time 
for corresponding 
temperature inte- 
gral (from case C) 


(N-N„)2 


Maix. temperature ■=■ 


m _T. 
a(Ng-No) 


Time constant - 


Rise time ■ 


Bjc. temperature 


integral 


i'lmo constant 


\ I !• I I I 
^^Cemperature integral 


?igure 4. - Various control parameters as functions of E for speed control with temperature 
limiting crlteriona (case B), when accelerating to limiting temperature. 
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Flgiire 7 . - Schematic diagram of control for turbojet engine for 
minimum rise time with temperatiore-llmlting criterions (case C) . 
Stability device must be added to this system. > 





Values of control parametesra 



• 6 1.0 i.2 1.4 1.6 1,8 2,0 


aJ^ (L-aNQ) 

Figure 8. - Various control parameters as function of for minimum rise time with 
temperature-limiting crlterlons (case C), when accelerating to limiting temperature. 
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